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Abstract 

In the first production run of the WASA experiment at COSY, the eta decay into three neutral pions was measured in proton- 
proton interactions at a proton beam kinetic energy of 1.4 GeV. The Dalitz plot of the three pions was studied using 1.2x10^ fully 
reconstructed events, and the quadratic slope parameter a was determined to be -0.027 ± 0.008(stat) ± 0.005(syst). The result is 
consistent with previous measurements and further corroborates the importance of pion-pion final state interactions. 
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1. Introduction 



The ?7 meson plays a special role in understanding low- 
energy Quantum Chromo Dynamics (QCD). Chiral symmetry, 
its realization in hadron physics at low energies and the role of 
explicit chiral symmetry breaking due to the masses of the light 
quarks (m, d, s) can be investigated using rf decays. 

The rj meson decays into three pions (77 — > ;r'^7r"7r" and 
Tj n^n^jfi) are among the major decay modes (the branching 
ratios are 32.6% and 22.7%, respectively IH) despite the fact 
that isospin is not conserved in these processes. It appears that 
these decays are driven by an isospin breaking term in the QCD 
Lagrangian which is proportional to the quark mass difference 
wid - Any contribution from electromagnetic processes 

is suppressed IJ, |4l |5|]. The decay is closely related to pion- 
pion scattering, an elementary low-energy QCD process. The 
lowest order contribution to the decay mechanism is given by 
Current Algebra (C A) 13, 01 ■ The predicted partial decay width 
is more than four times lower than the measured value. On the 
other hand, the distributions of the decay products are described 
within a few percent accuracy. 

The decay amplitude of tj into three pions can be described 
using the following two Dalitz variables: 
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1. 



(1) 



Here T, are the kinetic energies of the pions in the rest frame 
of the 77 meson, 3(r) = T\ + T2 + - m,, - iniyj where m„ is 
the pion mass and the tt", ji^ mass difference is neglected. The 
boundaries of the Dalitz plot are shown in Fig.[T] 

The decay amplitude predicted by CA is a linear function of 
y for the tj — » n^n^if' decay what implies that the rj — » ifn^Tf 
decay amplitude (^) is constant. Experimentally, a small devia- 
tion from a uniform Dalitz Plot density distribution is observed. 
The lowest term in the expansion about the center of the Dalitz 
plot is given by: 



|J^(z,0)r = c(l+2az) 
where z and (p variables are related to x and y via: 
X— Vicos0; y— -\fz sin <p. 



(2) 



(3) 



c is a constant factor and a is the quadratic slope parameter. 

The slope a was measured to be negative and small, which 
leads to a decrease of the Dalitz plot density at the border by a 
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Figure 1: Symmetrized Dalitz plot for the — > rfin^n^ decay. Dashed lines 
represent threshold for tt^tt^ — > tt^it* rescattering. 



few percent. The explanation of this effect poses a challenge for 
Chiral Perturbation Theory (ChPT), the effective field theory of 
QCD in the low-energy region. The leading order calculations 
in ChPT coincide with CA whereas the next-to-leading order 
(NLO) calculations |^ significantly improve the agreement for 
the partial decay width but predict a small positive value for 
a. The large uncertainty of the recently carried out next-to- 
next-to-leading order (NNLO) calculations does not allow 
to decide the sign of the slope. The origin of the non-uniform 
density distribution in the Si/r" Dalitz plot are in part pion-pion 
interactions in the final state. In contrast to the perturbative cal- 
culations, as soon as the rescattering effects are considered to 
infinite orders via unitarization of the decay amplitude using 
dispersion relations [10, 11] or iteration of the Bethe-Salpeter 
equation 11211 the sign of a turns out to be in accordance with 
experiment. The calculations indeed predict negative values for 
a in agreement with the experimental results which are sum- 
marized in Table [T] together with the different theoretical pre- 
dictions. 

Furthermore, 77 — > Stt decays are considered to be a source 
of precise constraints for the light quark mass ratios il3n . In 
a recent approach the constraints are derived entirely from the 
experimental partial decay width and from the data on 77 ^ 
7T°n°n° Dalitz plot slope ifljl . 

The first precise experimental determination of the a param- 
eter was carried out by Crystal Ball at AGS [15] in 2001, using 
a pion beam impinging on a liquid hydrogen target. In 2005, 
KLOE released their first result of a hig h statistics measure- 
ment of a (using — » O — > 777) 11611 which differed from 
the Crystal Ball result by three standard deviations. The CEL- 
SIUS/WASA IitIi experiment has confirmed the negative sign 
of the slope parameter a. However, the achieved accuracy did 
not allow to resolve the discrepancy between the Crystal Ball 
and the early KLOE results. Very recently, the KLOE data were 
reanalyzed and a new a value was presented ifisll that is consis- 
tent with Crystal Ball. However, the final results are not pub- 
lished yet. The Crystal Ball collaboration has recently collected 
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additional statistics at the MAMI facility that are currently be- 
ing analyzed lll9ll . 

The phase space of the rj — > Stt" decay covers the thresh- 
old for the t/^tt^ — » tt^tt' rescattering process as a consequence 
of the fact that the neutral pion mass is slightly lower than the 
charged pion mass. The boundaries correspond to the tt^tt'^ in- 
variant mass being equal to 2m;r± . In 2004 a cusp like structure 
in the tt'V invariant mass for the — > n^n^n^ decay was 
observed by the NA48/2 collaboration i22|]. The effect was pre- 
dicted already in 1997 by MeiBner et al. 12311 and the interpre- 
tation of the NA48/2 results was given by Cabibbo 0]. The 
process provides a precise determination of a combination of 
the 7 = and 7 = 2 pion-pion s-wave scattering lengths, aq and 
02- A similar effect was observed in the Ki — > 37r" decay 
and it should be present also in the rj — > 3;t" and rj' jf'n^rj 
decays. There are important consequences of this phenomenon 
for the analysis of the rj — > iifi decay [26^ ,27, 28, 17, 29]. The 
amplitude is no longer a function only of the z variable but has 
to explicitly depend on the Mandelstam variables. The descrip- 
tion of the amplitude as a polynomial function will fail in the 
cusp region. For example the dependence of the (p averaged 
ampUtude squared (|^P): 



1 C^" 

27r Jo 



l^(z,0)r# 



(4) 



on the z variable cannot be linear for 0.6 < z < 0.9. However 
the most sensitive variables to search for the cusp effect are the 
invariant masses of the pion pairs (or kinetic energies of the 
pions). They correspond to the projections onto the dotted lines 
in Fig.[T] 

The results presented in this paper are based on a measure- 
ment using proton-proton interactions with the WASA detec- 
tion system recently installed at the Cooler Synchrotron COSY 
at the Forschungszentrum Jiilich GmbH. 

2. The experiment 

The WASA detection system Q was transported to COSY 
(COoler SYnchrotron) in the Summer 2005 111. The COSY 
facility at Jiilich offers polarized and phase space cooled proton 
and deuteron beams with momenta up to 3.7 GeV/c l32l 13311 . 
The available beam momentum range allows the production of 
n, K, rj, <jj, and 77/ mesons well above the production thresholds 
in proton-proton and proton-deuteron interactions. Most of the 
final states in pN, pd, and dd reactions can be detected due to 
the nearly 4-n acceptance of the WASA detector for charged and 
neutral particles. 

The WASA detector |'30| was designed and optimized for 
studies of production and decays of light mesons in hadronic 
interactions. It consists of three main components: 
the Forward Detector - used for tagging and triggering on me- 
son production, the Central Detector - used for measuring neu- 
tral and charged meson decay products, and the unique pellet 
target system. The target beam consists of 30 fim diameter pel- 
lets of hydrogen or deuterium, providing a high target density 
in the order of 10'^ atoms/cm^. 




Central Detector 



FRH 

Forward Detector 



Figure 2: Schematic side view of tlie WASA detector setup at COSY. 



The Central Detector surrounds the interaction region and 
is designed for detection and identification of photons, elec- 
trons, and charged pions. It consists of an inner drift cham- 
ber (MDC), a superconducting solenoid providing the magnetic 
field for momentum measurements, a barrel of thin plastic scin- 
tillators for particle identification and triggering, and an elec- 
tromagnetic calorimeter The amount of structural material is 
kept to a minimum to reduce disturbance of the particles. The 
beryllium beam pipe wall is 1 .2 mm thick and the material of the 
superconducting solenoid corresponds to 0. 1 8 radiation lengths. 

The calorimeter (SEC) is used for detection and recon- 
struction of particles in the polar angle range from 20° to 
169°, which is about 96% of the geometrical acceptance. The 
calorimeter consists of 1012 sodium doped Cs I crystals. The 
trapezoidally shaped crystals with lengths from 20 to 30 cm 
correspond to ~16 radiation lengths. The energy resolution of 

5%/ V£[GeV] 



the calorimeter is o-{E)/E ^ 5%/ V^TGeV] OOll . 

The forward detector is designed for detection of particles 
emitted in the polar angles from 3° to 18°. It allows for identifi- 
cation and reconstruction of protons from the pp pprj reac- 
tion close to threshold. The precise track coordinates are given 
by four sets of straw proportional chambers (FPC). Kinetic en- 
ergies are reconstructed using the AE information in the layers 
of plastic scintillators of different thickness. In addition, the sig- 
nals are used for triggering. The kinetic energy of the protons 
can be reconstructed with a resolution of cr{T)/T ~ 1.5 - 3% 
for kinetic energies below 400 MeV. 

The detector setup is nearly the same as used in the previ- 
ous CELSIUS/WAS A experiment UtIi . The main change is a 
completely new readout system with charge-to-digital convert- 
ers based on a flash ADC l34ll35ll . In addition, the fourth layer 
of the forward range hodoscope (FRH) was replaced by two 
new thicker layers (15 cm instead of 1 1 cm). 

The rj mesons have been produced in proton-proton in- 
teractions at a beam proton energy of 1.4 GeV (momentum 
2.141 GeV/c). The beam energy corresponds to a center of 
mass excess energy of 56 MeV for the pp pprj reaction and 
the cross section is 10 fih. The results presented here are based 
on about 2.4 pb^^ of data collected during the first WASA-at- 
COSY production run and correspond to 94 hours of data tak- 
ing. At the trigger level, events with two tracks in the forward 
detector and at least two hit groups in the calorimeter with en- 
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Slope parameter 
a + Stat + syst 



Comment 



year Reference 



-0.022 ± 0.023 


GAMS-2000 


(1984) 


[20j 


-0.052 + 0.017 ±0.010 


Crystal Barrel 


(1998) 




-0.031 + 0.004 ±0.004 


Crystal Ball 


(2001) 


U5J 


-0.014 ±0.004 ±0.005 


KLOE, preliminary 


(2005) 




-0.027 ± 0.004 +0 0°^ 


KLOE, reanalysis 


(2007) 




-0.026 ±0.010 ±0.010 


CELSIUSAVASA 


(2007) 


uu 


-0.027 ± 0.008 ± 0.005 


WASA-at-COSY 


(2008) 


this result 





CA 


(1967) 




+0.015 


ChPT NLO 


(1984) 




-(0.007 + 0.014) 


ChPT NLO+dispersive 


(1995) 




-0.007 


UChPT 


(2003) 




-0.031 ±0.003 


UChPT fit 


(2005) 


{12} 


+0.013 + 0.032 


ChPT NNLO 


(2007) 


m 



Table 1 : Overview of experimental and theoretical results for the slope parameter a. 
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MM(pp) [GeV/c^ IM(ey) [GeV/c^ 

Figure 3: Missing mass of two protons, MM(pp) (left) and invariant mass of 
six photons, IM{2y) (right) after selection of the pp(yy final state. 

ergy deposits of more than 50 MeV were accepted. In addition, 
a veto on signals in the plastic barrel (PSB) was required, aim- 
ing to select only neutral particles. 

3. Data Analysis 

In the offline analysis, the pp6y final state is selected by re- 
quiring six hit clusters with energy deposit at least 20 MeV in 
the calorimeter consistent with photons and two proton tracks 
in the forward detector. The conditions provide a clean data 
sample of about 8-10^ events. In order to select candidates 



from the six reconstructed photons, all fifteen possible combi- 
nations of the three photon pairs were considered. For each 
combination, the quantity is calculated: 

>^j = h ^2 ' ; = 1.2, 3, ...15 (5) 

i 

where IM{i, i) is the invariant mass of the i-th yy pair for the 
j-th combination, cr - \\ MeV/c- is the resolution of the yy 
invariant mass. The combination with minimum value ofx^j is 
selected. Only events with;\^J < 15 and a missing mass of two 
protons between 0.535 GtWjc^ and 0.565 GeV/c^ are kept. Fig. 
[3] shows the missing mass of two protons (left panel) and the 
invariant mass of six photons (right panel) after applying the 
cut on the missing mass and after fitting the individual pions. 




0.2 0.4 0.6 0.8 1 

P(X2) 

Figure 4: The distribution of the kinematic fit probabiUty with the — > Sw" 
decay hypothesis (Pix^)) for the experimental data (solid line), Monte Carlo 
simulation of the signal (dashed line), and background (dotted fine). The verti- 
cal line represents the selection for the final DaHtz plot analysis. 
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In order to reconstruct the Stt" Dalitz plot, a kinematic fit 
is applied. In ref. 11711 . the full final state including the two 
protons was considered in the fitting procedure. Here, only the 
77 — > Stt" decay system is fitted, based on the reconstructed pho- 
ton angles and momenta. This approach reduces systematic un- 
certainty due to proton reconstruction. Different reconstruction 
uncertainties (as a function of energy and angle) were obtained 
from a full GEANT Monte Carlo detector simulation, with res- 
olution parameters matched to reproduce experimental distribu- 
tions. The experimental distribution of the kinematic fit prob- 
ability (Pix^)) is compared in Fig. |4]to the Monte Carlo sim- 
ulation. The simulation includes background from direct three 
pion production, pp ppTi'jf'jTp, the cross section was ob- 
tained by interpolation of the results from the CELSIUS/WAS A 
experiment ll36ll . The relative amount of background is less than 
4% in the final event sample after P{x^) > 0.1 cut (vertical line 
in Fig. nil. 

4. Extraction of the slope parameter 




Figure 5: Extracted dependence of on the z variable. Tlie solid line is a 
c(l + 2a7) fit. The dashed line is a prediction of the cusp effect f37ll . 

Based on the analysis procedure, the efficiency corrected ra- 
dial density distribution |^|- is shown in Fig.|5] The efficiency 
correction is obtained by dividing the measured z distribution 
by the result of the Monte Carlo simulation, assuming a - Q. 
A linear fit to the data points is applied to extract the slope pa- 
rameter a, yielding a statistical uncertainty of cr{a) - 8 ■ lO"-'. 
The achieved resolution in z is cr(7)=0.055. 

The systematic uncertainty of the result is estimated by vary- 
ing one by one all parameters that are important in the analy- 
sis. Table 12] summarizes some of these studies. The dominant 
contribution comes from an uncertain fraction of the interac- 
tions with gas stemming from pellet evaporation or with pellets 
bouncing in the beam pipe. This causes a spread in the vertex 
position. 

The combinatoric purity of the selected event sample depends 
on the cut on x^j probability of the best combination and on the 
confidence level of the fit. A variation of these parameters has 



only a minor effect on the obtained slope result. 
Another source of systematic errors is imposed by the remain- 
ing background. This effect is under control by changing the 
width of the proton-proton missing mass cut and, hence, vary- 
ing the relative amount of background in the range from 4 % to 
12 %. The observed effect contributes with 0.002 to the slope 
parameter uncertainty. 

The kinematic fit performance relies on the precise understand- 
ing of the errors of reconstructed quantities. This effect was 
studied by replacing the polar angle dependent parametrization 
of reconstruction uncertainties with a much simpler polar angle 
independent description, again showing only a small effect on 
the result (see table |2]i. 



Source of systematic uncertainty 



RMS 



Combinatoric background 0.001 

Missing mass cut 0.002 

Background Stt" production 0.001 

Vertex position 0.004 

Confidence level cut 0.001 

Error parametrization for the kinematical fit 0.002 



Overall systematic uncertainty 



0.005 



Table 2: The main contributions to the systematic uncertainty as discussed in 
the text. The overall systematic uncertainty was calculated as the square root of 
the summed squares of all contributions. 

Taking into account the systematic studies, the final result for 
the extracted slope parameter a is 

a = -0.027 + Q.Om(stat) ± Om5(sysf) 

based on 120000 events. The result agrees within one standard 
deviation with the high statistics measurements of the slope 
parameter performed by the KLOE collaboration llSll and the 
Crystal Ball collaboration 11511 . It is also consistent with the 
previous measurement performed by the CELSIUS/WAS A col- 
laboration I17il . The shape of the z distribution including the 
cusp that emerges from a virtual n^ji^ intermediate state fol- 
lowed by a 7r^;7r" — > tt'V transition, can be calculated from the 
rj — > n^n - jfi amplitude and the n - n scattering lengths. The 
cusp effect leads to a broad local minimum in the radial den- 
sity of the Dalitz plot for 0.6 < z < 0.9 due to the contribution 
of the regions where the invariant mass of the two pions less 
than Inij,* (see Fig.[T]i. The dashed line in Fig. |5] was obtained 
using a parameterization of the n — > ji^n^jf' decay amplitude 
from the KLOE collaboration 13811 and a nonrelativistic effec- 
tive field theory |2^ 37]. It is seen that the statistical precision 
of the data is insufficient to investigate the cusp effect. 



5. Outlook 

One motivation for the presented study was the striking dis- 
crepancy between KLOE and Crystal Ball results for the 77 
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3;r" Dalitz plot slope parameter a. Now, after the recent reanal- 
ysis of the KLOE data, the three major experiments focusing on 
the measurement of rj decays, KLOE at DAPHNE in Frascati, 
Crystal Ball now installed at MAMI in Mainz, and WASA-at- 
COSY in Jiilich, obtain consistent results. The non-zero value 
of the slope parameter is clear indication for the importance of 
final state pion-pion interactions. The theoretical understand- 
ing of these processes has advanced significantly, also triggered 
by the firm experimental situation. The Dalitz plot slope param- 
eter provides a very sensitive test of the ChPT predictions. The 
ChPT NLO and NNLO calculations that assume m,fi - m„± in 
the loops, indicate a positive sign of the slope parameter The 
uncertainty is however large and the negative sign is not ex- 
cluded iH. 

Among the experiments studying t] decays, WASA is 
presently the only experiment focusing on hadronic production 
of eta mesons in pp or pd scattering. The experiment com- 
bines the high production cross section especially in the pp 
reaction with the capability to run at high luminosities up to 
lO-'^cm^^i"'. This will allow to collect high statistics data sam- 
ples of J] mesons (10^ decays and more) with WASA-at-COSY. 
High statistics data samples provide further increased accuracy 
and sensitivity in the study of small effects, like the cusp struc- 
ture in T] -* 37r°. Recently WASA-at-COSY has collected a 
large data sample of pd — > ^Hct] events with low background 
and an unbiased trigger requirement, enabling studies e.g. of 
J] — > T^Tfin^ and 77 — > -k^tCti^ decays simultaneously. 
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